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We have measured the thermal conductivity of a TICUCI3 single crystal in magnetic fields 
up to 14 T. ft has been found that the temperature dependence of the thermal conductivity 
exhibits a sharp peak at 4 K in zero field, which is suppressed by the application of magnetic 
fields up to 7 T. The peak is concluded to be attributable to the enhancement of the thermal 
conductivity due to phonons because of the formation of a spin-gap state. In high magnetic 
fields above 7 T, on the other hand, another sharp peak appears around 4 K and this is 
enhanced with increasing magnetic field. This peak is regarded as being attributable to the 
enhancement of the thermal conductivity due to magnons and/or phonons because of the 
drastic extension of the mean free path of magnons and/or phonons in the Bose-Einstein 
condensed state. 
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Recently, the thermal conductivity in quantum spin systems has attracted considerable 
interest because it is closely related to the spin state. For example, in a spin-gap state, both 
the phonon-magnon scattering rate and the magnon-magnon scattering rate markedly de- 
crease, so that both Kphonon (thermal conductivity due to phonons) and /i magIlon (thermal 
conductivity due to magnons) are expected to be enhanced. Actually, enhancement of At p honon 
has been observed in the spin-gap state of the spin-Peierls system CuGeOs 1 ) and also of the 
two-dimensional orthogonal spin-dimer system SrCu2(B03)2- 2,3 ^ As for the enhancement of 
K ma gnom on the other hand, it depends on the magnitude of the bandwidth of the excited 
state of magnons whether K mag non is enhanced or suppressed. That is, in systems with a large 
bandwidth such as the two-leg spin-ladder system Sri4Cu2404i, 4 ~ 8 ) the enhancement has been 
observed, while the suppression has been observed in systems with a small bandwidth such 
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as CuGeC^. 1 ) In the antiferromagnetically ordered state of the two-dimensional spin system 
CU3B2O6, meanwhile, a marked increase in Kmagnon has been observed at low temperatures be- 
low the Neel temperature T^. 9 ^ Furthermore, it has been found that the thermal conductivity 
is sensitive to changes of the spin state in magnetic fields, such as reduction of a spin-gap and 
a spin- flop transition. 1_3 > 9 ) Accordingly, it is considered that the thermal conductivity mea- 
surement is a potential probe detecting a change of the spin state as a function of temperature 
and/or magnetic field. 

In this paper, we have applied the thermal conductivity measurement to a quantum spin 
system TICuCbj, which exhibits changes of the spin state, depending on temperature and 
magnetic field. The crystal structure of this material is monoclinic and composed of planar 
dimers of CU2CI6, 10 ) as shown in the inset of Fig. 1. The dimers are stacked on top of one 
another and form infinite double chains parallel to the a-axis. These double chains are located 
at the corners and center of the unit cell in the 6c-plane and are separated by Tl + ions. As 
for the spin state of TICUCI3, the ground state is a spin-gap state with the excitation gap 
A = 7.7 K in zero field. 10 ~ 12 ) From the analysis of the dispersion relation obtained by the 
inelastic neutron scattering experiment, it has been found that the spin gap of TICUCI3 is 
due to the strong antiferromagnetic interaction J = 5.68 meV in the planar dimer of CU2CL3 
and that the neighboring dimers are coupled by the strong interdimer interactions along 
the double chain and in the (1 2) plane. 13, 14 -* Our motivation of the thermal conductivity 
measurement for TICUCI3 arises from the three-dimensional magnetic ordering induced by 
the application of magnetic fields higher than the gap field H g ~ 5.5 T. 12 ) Under a constant 
magnetic field above ~ 5.5 T, the temperature dependence of the magnetization exhibits a 
cusplike minimum at Tn, and the ordering temperature increases with increasing magnetic 
field. 11 ) From a theoretical point of view, both the observed temperature dependence of the 
magnetization and the field dependence of the ordering temperature have been qualitatively 
well described in terms of the Bose-Einstein condensation (BEC) of magnons. 15 ) Actually, this 
suggestion has been confirmed by the observation of ordering of the transverse component of 
spins in the elastic neutron scattering experiment. 16 ) Moreover, very recently, the Goldstone 
mode characteristic of the BEC of magnons, theoretically suggested by Matsumoto et al., 17 ^ 
has actually been observed in the inelastic neutron scattering experiment by Riiegg et al.. 18 ^ 
In the BEC state, the behavior of the thermal conductivity is an attractive issue, because 
very large values of K m agnon are expected from the analogy of the thermal conductivity in the 
superfluid (BEC) state of liquid 4 He. 19 ) Therefore, in this study, we measured the thermal 
conductivity of TICUCI3 along the [2 1] direction, where the interdimer coupling is the 
strongest, and investigated the temperature and magnetic-field dependences. As a result, we 
observed drastic enhancement of the thermal conductivity in the BEC state for the first time. 

Single crystals of TICUCI3 were grown by the vertical Bridgman method. The details are 
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described in the previous paper. 12 ) Thermal conductivity measurements were carried out by 
a conventional steady-state method. One side of a rectangular single crystal (7.0 || [2 1] x 
1.0 x 1.0 mm 3 ) was anchored on the copper heat sink with indium solder. A chip resistance 
of 1 kf2 was attached as a heater to the opposite side of the single crystal with GE7031 
varnish. The temperature difference across the crystal (0.01-1.0 K) was measured with two 
Cernox thermometers (LakeShore Cryotronics, Inc., CX-1050-SD). The system was proved to 
be satisfactory from the thermal conductivity measurement of a standard reference material, 
Austenitic Stainless Steel. Magnetic fields up to 14 T were applied parallel to the 6-axis. 20 -* 

Figure 1 shows the temperature dependence of the thermal conductivity along the [2 
1] direction in magnetic fields parallel to the 6-axis up to 14 T. It is found that the thermal 
conductivity in zero field increases with decreasing temperature at low temperatures below 40 
K and shows a sharp peak at 4 K. With the application of magnetic fields up to 7 T, the peak 
is suppressed. In high magnetic fields above 7 T, on the other hand, the thermal conductivity 
increases with decreasing temperature and exhibits a kink and another sharp peak around 4 
K. This unusual field dependence of the thermal conductivity is clearly seen in Fig. 2(a). 

First, we discuss the origin of the sharp peak observed around 4 K in low magnetic 
fields below 7 T. Taking into account the theoretical calculation by Matsumoto 21 - 1 that the 
velocity of the excited state of magnons, namely, triplet magnons, is negligibly small in the 
spin-gap state of TlCuCkj, the contribution of « m agnon is neglected. Hence, at a glance, the 
temperature dependence of the thermal conductivity in zero field appears to be regarded as 
a typical behavior of Kphonon which exhibits a peak at a low temperature as a consequence 
of the freezing out of Umklapp processes of phonons with decreasing temperature on the 
one hand and the decreasing number of phonons on the other. However, the case is not 
simple, because Umklapp processes are generally independent of magnetic field so that the 
suppression of the peak in magnetic fields cannot be explained simply. It is likely that the 
peak is related to the formation of the spin-gap state, because the increase in the thermal 
conductivity with decreasing temperature below 40 K is in good correspondence with the 
decrease in the magnetic susceptibility below 40 K due to the formation of the spin-gap 
state. 12) As reported for CuGeCV) and SrCu 2 (B0 3 ) 2 , 2 ' 3) in fact, 

^phonon is enhanced in the 

spin-gap state, owing to the reduction of the phonon-magnon scattering. The peak is regarded 
as a result of the formation of the spin gap on the one hand and the decreasing number of 
phonons on the other. In magnetic fields, on the other hand, the suppression of the peak 
is interpreted as being attributable to the enhancement of the phonon-magnon scattering 
because of the reduction of the spin gap. Actually, the magnetic-field dependence of the 
peak in magnetic fields up to 7 T is similar to those in the spin-gap state of CuGeO-3 1 ) and 
SrCu 2 (B0 3 ) 2 . 2 ' 3) 

To prove this interpretation, we perform a fitting analysis for the magnetic-field depen- 



3/11 



J. Phys. Soc. Jpn. 



Letter 




J I I I I I I I I I I I I I I I I I I I I I I L 



10 20 30 40 50 

T(K) 

Fig. 1. Temperature dependence of the thermal conductivity along the [2 1] direction in mag- 
netic fields parallel to the 6-axis. The inset shows the crystal structure of TICUCI3. The dotted 
parallelogram indicates the unit cell. 

dence of the thermal conductivity at 4.90 K shown in Fig. 2(a). Usually, K p honon is given by 
tphonon = |CphononWp h onon T phonon, where C phonon is the specific heat of phonons, v pho non the 
velocity of phonons and r p honon the relaxation time of phonons. The scattering rate of phonons 
r phonon i s described as r~^ onoa = X)i r i -1 > with i denoting different scattering processes such 
as the phonon-magnon, phonon-boundary, phonon-defect and Umklapp scatterings. Changes 
of Chorion and Rphonon by the application of magnetic field are usually as small as at most 
1%, even if softening of phonons occurs because of the change of the spin state. 1-3 ' 22 ' 23 ) 
Therefore, the magnetic-field dependences of C pn onon and fphonon are neglected. In this simple 
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Fig. 2. (a) Magnetic-field dependence of the thermal conductivity along the [2 1] direction at 4.90 
K. Magnetic fields were applied along the 6-axis. The solid line denotes the best-fit curve of 
l/{a + fccxp(— A ^ )}. (b) Magnetic-field dependence of the mean free path of magnons along the 
[2 1] direction at 4.90 K. Error bars are due to errors of C ma gnon- 25 ^ The solid line denotes a 
fitting curve proportional to \H — H g \ n with n = 1.0. 



analysis, only the magnetic-field dependence of r p honon due to the phonon-magnon scattering 
is taken into account. The magnetic-field dependence of T~£ onoa is assumed to be simply due 
to the number of triplet magnons, so that the magnetic-field dependence of ftphonon is given 
by Rphonon (H) ~ l/{o + feexp(— A ^ )}, where a and b are constants and A(H) is the spin 
gap changing linearly with magnetic field. 12,13 ' 17 ) As shown in Fig. 2(a), the best fit curve is 
well overlapped on the experimental data points in low magnetic fields below 8 T. Thus, the 
suppression of the thermal conductivity with increasing magnetic field below 8 T at 4.90 K is 
well described in terms of the increase in the phonon-magnon scattering rate. Accordingly, it 
is proved that the sharp peak observed around 4 K in low magnetic fields below 7 T is due to 
the enhancement of Kphonon in the spin-gap state and that the suppression of the peak in low 
magnetic fields is attributed to the increase in the phonon-magnon scattering rate because of 
the reduction of the spin gap. 

Next, we discuss the origin of another sharp peak observed around 4 K in high magnetic 
fields above 7 T. Here, we define T K as the temperature where the thermal conductivity 
exhibits a kink, as indicated by arrows in Fig. 3. It is found that T K depends on the magnetic 
field H. T K values are clearly determined as the local minimum points in the k(H) — re(0) vs 
T plot and plotted in the H — T phase diagram in Fig. 4 together with BEC temperatures 
estimated from the magnetization 12 ) and the elastic neutron scattering 16 ) measurements. It is 



5/11 



J. Phys. Soc. Jpn. Letter 

found that T K is in good agreement with the BEC temperature. Furthermore, the inflection 
field (8 T at 4.90 K) in the magnetic-field dependence of the thermal conductivity shown in 
Fig. 2(a) is located on the phase boundary between the BEC state and the normal state, and 
hereby is regarded as H g at 4.90 K. Therefore, it is concluded that the origin of the sharp 
peak around 4 K in high magnetic fields above 7 T is due to BEC. 

Finally, we discuss whether the sharp peak around 4 K in high magnetic fields above 7 
T is due to the enhancement of Kphonon or K ma g non - As mentioned above, K ma g non is expected 
to bring about the enhancement on the basis of the analogy of the thermal conductivity in 
the superfluid state of liquid 4 He being much larger than that in the normalfluid state. 19 ) On 
the basis of the two-fluid model, 24 ) it is understood that the mean free path of the remain- 
ing normalfluid carrying heat becomes drastically long in the superfluid state, because the 
normalfluid-superfluid scattering rate is very small. In the case of TlCuCi3, by analogy, the 
mean free path of the normalfluid magnons may be expected to become drastically long in 
the BEC state so that Av ma g non is drastically enhanced. At much lower temperatures below 
the BEC temperature, on the other hand, /t ma gnon should decrease with decreasing tempera- 
ture, owing to the decreasing number of normalfluid magnons. This analogy may be plausible, 
but it is not sure whether such a hydrodynamic heat conduction due to the mutual conver- 
sion of the normalfluid and superfluid in liquid 4 He is valid in TICUCI3 as well. However, 
it is worthwhile examining the extension of the mean free path of magnons, ^magnon? in the 
BEC state of magnons. In the following, we estimate the magnetic-field dependence of Zmagnon 
from the magnetic-field dependence of the thermal conductivity at 4.90 K shown in Fig. 2(a). 

^magnon(-f^) is expressed as Kmagnon(-fO — Cmagnon (-^O^magnon (-fO^magnon (-^0 > with Cmagnon(^) 

being the specific heat of magnons and w m agnon(-£0 the velocity of magnons. Assuming that 
the thermal conductivity at H g = 8 T is equal to K p honon and that At p honon is independent of 
magnetic field above H g , rc m agnon(-fO is approximately estimated to be k(H) — n(H s ). On the 
other hand, C ma g non (-£0 is approximately estimated to be C(H) — C(0) from the specific heat 
measurements by Oosawa et al., 25 ^ and v m3jgnon (H) has been theoretically calculated by Mat- 
sumoto. 21 ) The values of calculated l magn0 n(H) are thus shown in Fig. 2(b). Error bars are due 
to errors of C magn0 n(-ff), because C(0) includes the Shottky-type specific heat of magnons. 25 ) 
It is actually found that Zmagnon exhibits a rapid increase with increasing magnetic field above 
Hg. ^magnon increases in proportion to \H — H g \ n with n = 1.0 around H g and is as large as 
7.5 x 10 4 A at 14 T. Above H g , t> ma gnon also increases with increasing magnetic field, but the 
value at 14 T is at most 2.0 x 10 5 cm/s 21 ' which is not so different from the typical value 
of Uphonon- 22 ' 23 ^ On the other hand, the estimated value of Zmagnon at 14 T is about ten times 
larger than those of one-dimensional spin systems CuGeOs, 26 ) SrCu02 27 * ) and S^CuOs 27 ) 
where the large contribution of K mSLgnon has been pointed out. Such a large value of Zmagnon 
in TICUCI3 may be characteristic of the BEC state of magnons, suggesting the drastic sup- 
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pression of the magnon-magnon scattering rate as discussed on the basis of the analogy of 
the superfluid state of liquid 4 He. Accordingly, it is likely that the observed sharp peak of 
the thermal conductivity around 4 K in high magnetic fields above 7 T is due to the drastic 
enhancement of K magnon , owing to the increase in t> ma gnon and the drastic extension of /magnon- 

On the other hand, it is possible that the sharp peak around 4 K in high magnetic fields 
above 7 T is due to the enhancement of Kphonon- In this case, the origin of the peak is attributed 
to the extension of the mean free path of phonons in the BEC state where the phonon-magnon 
scattering is suppressed as the number of the normalfluid magnons decreases. One may claim 
that Kphonon in the BEC state cannot be larger than K pn0 non in zero field around 4 K, because 
the phonon-magnon scattering rate in the spin-gap state seems to be always smaller than 
that in the gapless state such as in the present BEC state. However, this may be misleading, 
because the phonon-magnon scattering in zero field may not be suppressed completely around 
the peak temperature (4 K) which is comparable with A = 7.7 K. To be more conclusive, 
thermal conductivity measurements for the Mg-substituted TlCuCbj are under way. 

In conclusion, for a TlCuChj single crystal, we have measured the thermal conductivity 
along the [2 1] direction in magnetic fields parallel to the 6-axis up to 14 T. It has been 
found that the thermal conductivity is drastically enhanced with increasing magnetic field in 
the BEC state. The enhancement has been regarded as being attributable to the enhancement 
of the thermal conductivity due to magnons and/or phonons because of the drastic extension 
of the mean free path of magnons and/or phonons. To the best of our knowledge, the present 
work is the first report on the enhancement of the thermal conductivity in the BEC state of 
magnons. 
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Fig. 3. Low-temperature part of the temperature dependence of the thermal conductivity shown in 
Fig. I. Each plot in a magnetic field is shifted vertically by 20 mW/Kcm for clarity. Arrows 
indicate the temperature below which another peak develops with decreasing temperature. 
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Fig. 4. Phase boundary between the Bose-Einstein condensed state and the normal state in TlCuCbj, 
determined from the thermal conductivity, magnetization 12 ) and neutron scattering 16 ) measure- 
ments. Magnetic fields were applied parallel to the 6-axis. The dashed line is a guide for the 
eyes. 
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